###########################CARREGAR PACOTES##############################
#########################################################################
require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(agricolae)
require(dae)

#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")



#########################################################################
#######################Planejamento do DIC###############################
#########################################################################
t=3 ##numero de tratamentos
r=4 ##Numero de repetições
n=r*t
DIC.parc=list(Repeticao=r, Parcela=t)
DIC.nest=list(Parcela="Repeticao")
tratamento=factor(rep(c(1:t), times=r))
DIC=fac.layout(unrandomized = DIC.parc,nested.factors=DIC.nest,randomized = tratamento , seed = 105)
layout=matrix(DIC[,5],t,r)
layout


#########################################################################
############################Análise do DIC###############################
#########################################################################


##Importar dados
dados=readXL("fungos.xlsx",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)


####GERAR MODELO
model=aov(vr~trat)
anova(model)
cv.model(model)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)


#Teste de Kolmogorov-Smirnov - NORMALIDADE
m=mean(e)
dp=sd(e)
ks.test(e,"pnorm",sd=dp)

#Teste de Levene - Homogeneidade de Variância
leveneTest(model)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model)